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ABSTRACT 

We examine the gas and stellar metallicities in a sample of HII galaxies from the Sloan 
Digital Sky Survey, which possibly contains the largest homogeneous sample of HII 
galaxy spectra to date. 

We eliminated all spectra with an insufficient signal-to-noise ratio, without strong 
emission lines, and without the [OH] A3727 A line, which is necessary for the deter- 
mination of the gas metallicity. This excludes galaxies with redshift ^ 0.033. Our final 
sample contains ^700 spectra of HII galaxies. 

Through emission line strength calibrations and a detailed stellar population anal- 
ysis employing evolutionary stellar synthesis methods, which we already used in pre- 
vious works, we determined the metallicities of both the gas and the stellar content of 
these galaxies. 

We find that in HII galaxies up to stellar masses of 5 • 10^ M©, enrichment mecha- 
nisms do not vary with galactic mass, being the same for low- and high-mass galaxies 
on average. They do seem to present a greater variety at the high-mass end, though, 
indicating a more complex assembly history for high-mass galaxies. In around 23 per 
cent of our HII galaxies we find a metallicity decrease over the last few Gyr. Our 
results favour galaxy evolution models featuring constantly infalling low-metallicity 
clouds that retain part of the galactic winds. Above 5 • 10^ Mq stellar mass, the 
retention of high metallicity gas by the galaxies' gravitational potential dominates. 

Key words: ISM: abundances - Galaxies: evolution - Galaxies: starburst. 
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1 INTRODUCTION 

HII galaxies are characterised by prominent emission lines. 
In fact, they are defined as having strong H/3 and QUI lines, 
but other ionised hydrogen (Hq, H7, etc.) and high excita- 
tion metallicity lines (Oil, Nil, and others) are strong too. 
These emission lines are driven by photoionisation of inter- 
stellar ga£_by_Jio^jimssive_st stellar popula- 
tions (|Sargent fc Searlelll970l : iFrenchlllQSdl '). The star form- 
ing activity in HII galaxies is so strong that it can certainly 
not have been m aintained at its present level durin g a Hub- 
ble time (see e.g. ISearle. Sargent fc BagnuoldllQT^ '). In fact, 
they are the most extreme case of star forming galaxies, 
showing the highest excitation emission lines, and thus the 
highest (relative) star formation rates (SFRs). 
They are also the lowest metallicity galaxies of the in- 
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ters tellar medium, with metallicities around one tenth so- 
lar (IVflchez fc Iglesias-Para,mdl2003l:|Perez-Montero fc Dfaz 



'200^; 'Kniazev et al. 2004": 'Kehrig et al.' '2006': 'Izoto v et al 
|2006,; Hagcleet aL 2008; P crcz-Montcro ct al. 201^, and 
others). 

It is now clear that the hypothesis of HII galaxies as 
galaxies experiencing their first star formation burst can 
be ruled out, and that their stella r content is predom- 
inantly ol d, older than one G yr dR aimann et al.l 
Cid Fernan dcs. Lca o fc Lacerdal 20031: iKong et al.l 
Westera ct al. 2004; 'Ci d Fernandes et al. I I2007I : lAsari et al.l 
2007; Hoyos ct al. 2007, and others). Nevertheless, both the 
high relative SFRs and the low metallicities indicate that 
HII galaxies are among the least evolved galaxies in exis- 
tence, which assi gns to them a special role as fossil record of 
galaxy evolution l|Leaueux et al.|[l979l ). They are, therefore, 
ideal objects for chemical evolution studies. 

Several authors have developed theoretical evolution- 
ary models of dwarf irregular galaxies, of which HII galax- 
ies are a sub-category. Some of these models predict the 
galaxies' present-day stellar - and gas metallicities, among 
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other properties. In these models, several chemical enrich- 
ment processes are i dentifi ed. 

iMacLow fc Ferraral (|l999l ) developed single-phased hydro- 
dynamical models for dwarf galaxies. They model the effects 
of repeated supernova type II (SNII) explosions from star- 
bursts on the interstellar medium, that is, the enrichment 
by the ejected metals on the one hand, and the gas loss 
by supernova winds on the other hand, taking into account 
the gravitational potential of their dark matter haloes. They 
find that, in galaxies with gas masses below 10^ Mq, most of 
the gas is blown away, and in galaxies with masses between 
10^ Mq and 10^ Mq, mainly the newly- formed metal-rich 
gas is blown away, whereas the already present, metal-poor, 
gas is retained. They do not calculate the final (present- 
day) metallicities of their model galaxies, but it is obvious 
that they remain low-metallicity, at least the ones with gas 
masses below 10^ Mr. 

iRecchi et al] (|200ll . l2002l ') simulate IZwl8, one of the lowest- 
metallicity galaxies known, using a single star formation 
burst model and a model with a doubly peaked star forma- 
tion history. They use a lower SNII heating efficiency than 
MacLow et al. in their models, but include type la super- 
novae (SNela), for which they use a higher heating efficiency. 
As a consequence, alpha elements are ejected less efficiently 
than in MacLow et al.'s model, resulting in higher metallic- 
ities. The iron-peak elements, which are produced in SNela, 
on the other hand, are still ejected, causing high [a/Fe] ra- 
tios. The (gas) metallicity in their single burst model re- 
mains very low, [O/H] S —1.3, depending on the galaxy 
mass. In their (more likely) double burst model, [O/H] can 
reach values up to -1, depending on the details of the mod- 
els, that is, the stellar yields used, the stellar initial mass 
function (IMF), the time between the bursts, the duration 
of the second burst, the gas density and - metallicity in the 
star forming region, and the total mass of the stars formed in 
the second burst. The final metallicities of the double burst 
model galajcies are between 0.6 and 1 dex higher than after 
the first burst. 

iTenorio-Tagle et al.l (|2006l ) include photoionisation and clus- 
ter wind in their two-dimensional hydrodynamic calcula- 
tions. They identify two mechanisms: The storage of clouds 
into a long-lasting ragged shell inhibiting the expansion of 
the thermalised wind, and the steady filtering of the shocked 
wind gas through channels carved within the cloud layer. 
They conclude that both processes must be at work in HII 
galaxies. Unfortunately, they make no prediction about the 
present-day metallicity of their model galaxies. 
In a series of papers, Recchi and collaborators use a new 
generation of chemo-dynamical m odels, some o f whic h in- 
clude inf ailin g cloud s. The first two. lRecchi et al.l l|2004l ') and 
I Recchi et al.l l|2006l ). are aimed at reproducing the proper- 
ties of two specific objects, IZ wl8 and NGC1569, respec- 
tively. The third one, .Rec chi fc HensleJ (HoOTi) 

investigates 

the infiuence of several factors, such as the cloudiness of the 
gas distribution, the IMF slope, and stellar yields. They find 
that models with continuous low level star formation peri- 
ods in the past followed by a quiescent phase and a recent, 
stronger burst best reproduce the chemical properties of the 
studied galaxies. They also find that, for a homogeneous gas 
distribution, metals get blown away by galactic winds, but 
cool infalling intergalactic clouds can hamper these galac- 



tic winds. The final metallicities in the their various models 
vary from -2.5 to -0.8. 

Although metallicity determinations of the various com- 
ponents of HII galaxies cannot discriminate, how much each 
of these mechanisms contribute to the chemical evolution, 
they can quantify the metallicity at various stages of a 
galaxy's evolution, i.e. at the time of the formation of the old 
(> 1 Gyr) stellar populations and at present, and thereby 
help to rule out some of the evolutionary scenarios, and sup- 
port others. HII galaxies have the advantage that it is possi- 
ble to determine independently the metallicities of the young 
and the old populations from their spectra. As the young 
populations were formed recently out of the same gas that 
is responsible for the emission lines, their metallicity can be 
assumed to be more or less the same as the one of the gas, 
which can be determined from the emission line strengths. 
The metallicity of the old populations, on the other hand, 
can be determined from the continuum and absorption fea- 
tures through population synthesis. 

The main goal of the present work is to measure the chem- 
ical enrichment in a homogeneous sample of HII galaxies, 
and examine possible trends with galactic properties such 
as total mass. This way, we can evaluate the different mod- 
els and make statements on the importance of the various 
enrichment mechanisms. 

Until recently, the discovery of HII galaxies was limited 
to the visual inspection of objective-prism surveys, intro- 
ducing ill-c ontrolled biases and selection effects (for compi- 
lation s, see lTerlevich et ah l ll99ll : lKehrig7 Telles fc Cuisinien 

|2004| . and references therein). 

The Sloan Digital Sky Survey (SDSS; lYork et al.l |200G| ) 
presents for the first time a c omprehensive databas e of galac- 
tic spectra. Data release 7 (jAbazaiian et al.ll2009l ) contains 
over 900'000 galaxy spectra that have been selected on clear 
magnitude limiting criteria. In an unprecedented manner, 
the SDSS allows to study star forming emission line galax- 
ies based on clear and quantitative criteria, and not mere 
visual inspection. 

In this work, we define clear emission line strength cri- 
teria to distinguish galaxies containing strongly excited gas 
from other emission line galaxies. We only include spectra 
that contain the [OH] A3727 A emission line, as this line is in- 
dispensable for an adequate determination of the gas metal- 
licity. This criterion excludes galaxies with redshift S 0.033. 
This way, we selected a homogeneous sample of ~ 700 HII 
galaxy spectra of high quality from the SDSS. 
In this sample, we determine the metallicities of both the 
gas and the old stellar content using emission line calibra- 
tion and stellar population synthesis methods, respectively. 
It is the first time an independent metallicity determination 
of the gas and the stars is performed on a clearly defined 
HII galaxy sample. 

The layout of this article is the following: In Sect. [S] 
the catalogue of spectra analysed in this work is presented. 
Sect.|3]gives a detailed description of the method we used to 
analyse the spectra. Subsection 13. II describes, how we deter- 
mined the gas metallicity, and Subsection 13.21 is dedicated 
to the determination of the stellar metallicities. The main 
results are given in Sect.|4]and discussed in Sect. [5] A sum- 
mary and the main conclusions can be found in Sect. [6] 
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Figure 1. Illustration of the emission line strength measurement 
(here H/3 in 51957-0273-418). In the left panel, the solid line 
represents the observed spectrum (rebinned to the BC03 wave- 
length grid), whereas the short-dashed line shows the best fit us- 
ing the BC03 high resolution library. The shaded region between 
these two lines shows the area used to calculate the emission-line 
strength. The right panel shows the emission line after subtract- 
ing the best fit. The shaded area corresponds to the shaded area 
between spectrum and best fit from the left panel. 



2 THE DATA BASE 

For this work, we used spectra from the SDSS data release 
7, which have a signal-to-noise of S/N > 4 per pixel at 
nig = 20.2, but usually much higher, and cover a wavelength 
range of 3800 to 9200 A. Most important of all, it contains 
a large number of galaxy spectra, that is around QOO'OOO. 

Any analysis however relies on previous clear crite- 
ria to define what HII galaxies are. First, we rejected all 
spectra whose wavelength ranges do not include the region 
from 3677 A to 3775 A, which is needed to measure the 
[Oil] A3727 A emission line strength, used in many strong 
line methods to determine the gas metallicity. This excludes 
galaxies with redshifts below 0.033. Another prerequisite for 
measuring this line strength is a sufficient signal-to-noise 
ratio in this wavelength region. We only kept spectra with 
5/A''[oii]3727 ^ 2.8. This criterion guarantees sufficiently low 
noise for the entire spectral range used in this work, since at 
longer wavelengths, the signal-to-noise ratio of the spectra 
is higher than around the [Oil] A3727 A line, i.e. Z 8 in g. 
The spectra were then corrected for Galactic (foreground) 
gas extinction using the value s given in the SDSS database, 
which were derived from the Schlegel. F inkbeiner fc David 
(|l998h reddening maps. Usin g those extinction values and 
the Galactic extinction law of iFitzpatrickl l|l999l ). we dered- 
dened the spectra. After that, we deredshifted them, and 
measured the emission line strengths. 

In order to properly determine the emission line strengths, 
we first had to remove the contribution from the absorp- 
tion lines of the underlying stellar populations. This was 
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Figure 2. Excitation diagram for all galaxies of our sample (black 
crosses). The grey crosses represent the test sample, whereas 
galaxies identified as AGNs are sh o wn a s dashes. The solid line 
shows models from iKewlev et al. 1 ll200ll ') used to separate the 
AGNs from star forming galaxies. 



done by subtracting high resolution spectra representing 
the stellar continuum from the empirical spectra as illus- 
trated in Fig [1] These high resolution spectra were com- 
piled by fitting the spectra of composite stellar populations, 
made up of three single stellar populations (SSPs), a young 
(^10 Myr), an intermediate-age (20 - 500 Myr), and an old 
one (5 Gyr), from the "BC03" integrated spectral energy dis- 
tribution (ISED) library to the continua of the galaxy spec- 
tra. The "BC03" library was produced using the Bruzual 
and Chariot (2003) Galaxy Isochro ne Spectral Synthesis 
Evolution Library (GISSEL) code (ICharlot fc Bruzual A.I 
Il99ll : iBruzual A. fc Charlotlll993l, 12003 ) implementing the 
Pado va 1995 isochrones ( Fagotto et al .' 1994: Girardi et ajj 
Il996l ) combined with the STELIB ((l7e Borgne et al., ,200^) 
stellar library. Since it is made up of empirical high resolu- 
tion spectra, the "BC03" library reproduces absorption line 
shapes well, and is ideal for this purpose. Table [1] lists all 
the lines we measured. 

However, incomplete coverage of the stellar parameter space 
and calibration uncertainties of the STELIB library result 
in systematic errors in the overall spectral shapes. Therefore 
the "BC03" library is not a good choice for determining the 
properties of the stellar content from a full spectral fit. For 
the stellar population synthesis, we perform another fit us- 
ing a different library, described in the next section. 
Finally, w e corr ected the spectra for internal gas extinction 
using the Calzetti (2001) attenuation law. The extinction 
constants E{B ~ V) = 0.44 x 6.60 log(7Hc«//H/3/2.87)/4.04 
were estimated from the Ha/H/3 Balmer decrements foUow- 
i ngX^ alzctti (200lj), adopting intrinsic ratios Ina/Inp = 2.87 
(iOsterbrocklll989l ). 

The factor of 0.44 stems from the correction for differen- 
tial extinction between the stellar populations and the gas 
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Table 1. Emission lines measured in this work. 



line 


central wavelength 


[Oil] A3727 A 


3727 A 


He 


3970 A 


H<5 


4102 A 


H7 


4340 A 


[OIII] A4363 A 


4363 A 


[Hell] A4686 A 


4686 A 


H/3 


4861 A 


[OIII] A4959 A 


4959 A 


OUT A5nn7 A 


5007 A 


[Hel] A5876 A 


5876 A 


[OI] A6300 A 


6300 A 


[SIII] A6312 A 


6312 A 


[Nil] A6548 A 


6548 A 


Ha 


6563 A 


[Nil] A6584 A 


6584 A 


[SII] A6717 A 


6717 A 


[SII] A6731 A 


6731 A 


[Oil] A7319 A 


7319 A 


[Oil] A7330 A 


7330 A 



following ICalzetti. Kinney fc Storchi-BergmannI l| 19941 ) . As 
our sample consists of HII galaxies only, which is a subclass 
of starburst galaxies, using this factor is justified. The emis- 
sion line strengths were dereddened as well, using the same 
E{B — V) values but without multiplying them by 0.44. 

To make sure our sample contains only HII galaxies, we 
limited it to galaxies with strong emission lines at a high 
excitation level, by applying the two criteria: 



log 



/ [OIII] A5007 A 



H/3 



^ 0.2 II log 



NII]\ 



Hq 



-0.5, 



(1) 



whereas [Nil] = [Nil] A6548 A + [Nil] A6584 A. Of course, 
these criteria are more restrictive than the traditional def- 
inition of HII galaxies as galaxies showing strong Balmer 
lines. In fact, as our criteria are designed to identify galax- 
ies with gas at a high excitation level, our sample is biaised 
towards extreme case HII galaxies. As a consequence, we 
expect the galaxies of our sample to show typical proper- 
ties of HII galaxies, but at a more pronounced level. For 
example, it consists of HII galaxies with low mass-to-light 
ratios and thus contains mainly dwarf galaxies. The reason 
for this approach, apart from emphasising on typical HII 
galaxy properties, is, to make sure that no "contaminat- 
ing" objects, such as LINERs, enter the sample. Then, we 
separated the star forming galaxies from AGNs, using the 
iKewlev et al.l (|200ll ) excitation line criterion. 



, / [OIII] A5007 A 
^'^n 



aei 

log([NII]/HQ) -0.47 



+ 1.19, (2) 



as shown in Fig. [2] 

Our final sample contains 712 HII galaxy spectra. The rela- 
tively low number of galaxies compared to the total number 
of over 900'000 galaxies in the SDSS data release 7 is due 
to our high requirements towards the quality of the spectra, 
and to the criterion that the spectra must contain all the 
lines necessary for our study. This way, we guarantee that 
the sample is not contaminated, and that the obtained re- 



sults are trustworthy. However, the sample is large enough 
to be statistically significant. 

We also defined a test sample, containing galaxies with 
some line emission, but not necessarily at a high level, that 
is Ha, H/3, [OH] A3727 A, [OIII] A4959 A, [OIII] A5007 A, 
[Nil] > 0. The reason for defining a test sample is, on the one 
hand, to identify properties typical of HII galaxies through 
comparison with the test sample, and, on the other hand, 
to check our method, by comparing the results for our test 
sample with the ones obtained in other studies using similar 
galaxy samples. As we needed to measure the [Oil] A3727 A 
line strength for the test sample as well, it is also limited 
to galaxies with redshift ^ 0.033. To guarantee sufficient 
quality, we applied the signal-to-noise criterion S/Ng ^ 5, 
where the S/Ng values were taken from the S DSS database. 
Here too, we removed the AGNs using the iKewlev et al.l 
l|200lh criterion. This way, our test sample contains 74989 
spectra. Given the criterion that the test sample galaxies 
show some line emission, it consists of star forming galaxies, 
probably mostly spirals. 

Fig. [2] shows both our sample of HII galaxies (in 
black), and the test sample (in grey) in a Baldwin-Phillips- 
Terlevich (BPT, [Baldwin et al.l Il98ll ) classification dia- 
gramme. 



3 METHOD 

3.1 Determination of tlie gas metallicities 

We first measured the gas metallicity, which corresponds to 
the metallicity of the young stellar populations. In order to 
derive gas abundances from its line emission, hydrogen lines 
are needed, and the lines of at least one ion, generally oxy- 
gen, in its two dominant ionisation stages. Unfortunately, 
the electron temperature method using the [OIII] A4363 A 
emission line could not be used as this line and other auroral 
lines are much too weak to be measured with the necessary 
precision at the signal-to-noise ratios of the SDSS spectra. 
In most spectra, we do not even detect the [OIII] A4363 A 
line at all. Hence, we are restricted to strong line methods. 
The [OH] A7330 A line could be used, but it is too sensitive 
to the electron temperature and is possibly "contaminated" 
by recombination contributions. 

We compare the gas metallicities derived from 11 different 
strong line methods. 

A frequently used i ndicator is the [Nil] A6548 A or the 
[Nil] A6584 A hue (|Denicol6 et al.1 12002| : iPettini fc Pagell 
I2OO4I ). However, since nitrogen is produced and destroyed 
during both primary and secondary nucleosynthesis, its 
abundance correlates in a non-evident way on (oxygen) 
metallicity. Therefore, it should be used only for rough es- 
timates, to discriminate between the different branches of a 
multiple-valued method, for instance. 

One su ch multip l e-valu ed indicator is the R23 param- 
eter (iPagel et al.1 1 19791 ). which sums up the fluxes of 
the [OH] A3727 A and two strong [OIII] lines, at 
4959 A and 5007 A. As -R23 is aproximately propor- 
tional to the oxygen abundance at low metallicities, but 
decreases at high metallicities due to cooling, this in- 
dicator has a low - and a high metallicity branch. In 
spite of this double- valuedness, the R23 parameter is 
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Figure 3. Comparison between the metaUicities determined using thelPii vugin &: ThuarJ ll2005h metho d (x-axes) and the ones deter- 
mined by nine other strong-line methods (y-axes). K. & K. 2004 stands for lKobuInickv Sz Kewlevi ||2004| ) . The dashed lines show unity, 
[0/H]x = [0/H]p For further details, see to the text. 



used in many met a llicity ca l ibratio ns (iMcGaughl Il99ll: 
Zaritskv et all Il994l: rPilvuginl l2001allbl: IXe^lev fc Dopital 



20021 : iKobulnicky fc Kewle vil2004l: iPilvurin fc Thuanll2005H . 

some of which include a recipe, how to determine, on which 
branch a given galaxy lies, whe reas others leave it up to the 
user to decide this. In one case. IZaritskv et al.l l| 19941 ) . a cal- 
ibration is determined only for the high metallicity branch. 
Another point that has to be taken in to account (although 
in some methods this is ignored, i.e. IZaritskv ctal] Il994l ') 
is, that R23 depends not only on the metallicity of the 



ionised gas, but also on the hardness of the ionising radi- 
ation, which can be quantified by the ionisation parameter 
q, representing the ionising photon flux pe r unit area divided 
by the number density of hydrogen a toms. iKewlev fc Dopit"al 
(.2003 ): lKobulnickv fc Kewlevi (|2004 ' 'new paramatrisation" 
and "best estimate" methos) use the [OIII]/[OII] ratio to cal- 
culate q in an iterative way from theory, and then determine 
[O/H] from R23 and q, whereas iMcGaughl l|l99ll ) uses the 
[OIII]/[OII] ratio dir ectly in an empirical calibratio n with- 
out determining q. In lPilvugirn (,2001a .b): iPflvugin fc ThuanI 
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l|2005ll , the so-called excitation or ionisation or, simply, P pa- 
rameter is introduced for their empirical calibrations. It is 
defined as the [OIII]/([OII]-f [OIII]) ratio, so the [OIII]/[OII] 
ratio is related to it, [OIII]/[OII] = P/(l - P). Recapitulat- 
ing, in all of these methods the gas metallicity is calculated 

from the R2Z parameter and the [OI II] /[Oil] ratio. 

Other strong line methods include iKobulnickv fc KewlevI 
l|2004 [Nil] method), com bining the [Nil] A6584 A line and 
the [OIII]/[OII] ratio, and iPettini fc Pagell l|2004l 'l. employ- 
ing the ratio between the [OIII] A5007 A and [Nil] A6584 A 
lines. 

Fig. [3] shows the metallicities of our HII galaxy sam- 
ple as determined by nine of the methods we studied, 
plotted against t h e one s using by the tenth method, 
IPilvugin fc ThuanI l|2005h . which is the one we ended up 
adopting. The on ly calibration no t show n here due to space 
limitations is the lPettini fc Page! Hooi) [Nil] one, which is 
a vertically squeezed and shifted downward version of the 
iDenicolo et al.l (|2002l ) calibration (first panel). 
Most methods yield similar relative metallicities, that is, 
the methods agree upon the metallicity differences between 
galaxies. However, when it comes to absolute [O/H] val- 
ues, systematic differences of up to 0.5 dex betwe en meth- 
ods ca n b e observed. All methods apart from IPilvugii] 
( 2001a, b); Pilvugin fc ThuanI l|2005l . x axis) result in un- 
realisticall y high gas metall i cities f or HII galaxies. Espe- 
cially the iKewlev fc Dopital l|2002l 'l: IKobulnickv fc KewlevI 
l|2004 ) ones hardly ever give metallicities belo w -0.5. On top 
of tha t, the lKewlev fc D opita ( 2002); Kobulni ckv fc KewlevI 
l|2004l . [Nil] and "new paramatrisation" ) methods did not al- 
ways c onverge. The horizontal lines in the lKewlev fc Dopital 
(|2002l ) panel are due to some cases, in which the expression 

Iog(0///) + 12 = gave complex val- 

ues. In these cases, we used only the real part, —ki/2k2, 
which results in identical metallicity values for galaxies 
within the same ionisation parameter range. When ignor- 
ing these galaxies, we obtain a similar relation as for the 



Kobulnickv fc KewlevI ||2004| ) ' 'new paramatrisation of the 



Kewlev fc Dopital (|2002l ')" method 



The method that yields the most realis tic metallicities for 
the ga laxies of our sample is the one bv lPilvugin fc ThuanI 
1(2003). It also accounts for the hardness of the ionising 
radiation by means of the P parameter. Therefore, we 
adopted it for our gas metallicity determination. To de- 
termine, on which branch a given galaxy lies, we used the 
me tallicity estimate based on the [Nil] A 6548 A/Ha ratio 
by IDenicolo. Terlevich fc TerlevichI ( 200d). as suggested by 



vario us authors ( Perez-Montero fc Diaj l2005l : IHovos et al.l 
I2OO7I . and others). Most galaxies of our sample lie on the 
high-metallicity branch. 

iLopez-Sanchez fc EstebanI (|2010l ) come to the same conclu- 
sion after comparing a series of empirical methods, like the 
ones presented in this section, with the results of the elec- 
tron temperature method for a sample of 31 Wolf-Rayet 
galaxies. They conclude, that the nowadays best suitable 
method for star-forming galaxies, where auroral lines such 
as JOIII] A4363 A are not observed, is IPilvugin fc ThuanI 
(|2005l ). and that other methods based on photionisation 
models yield metallicities system atically 0.2-0.3 dex higher 
and higher dispersion than the IPilvugin fc ThuanI (|2005l ') 
calibration. 



3.2 Determination of the stellar metallicities 

Then, we determined the metallicity of the stellar con- 
tent, that is, of the stellar populations making up the 
galajcies, using a pop ulati on synthesis method , described in 
ICuisinier et al.l (120061 ) and lLisker et al.l (120061 ). based on the 
full spectral fitting in the 3820-8570 A range whereas the 
upper limit can be lower than 8570 A (in some rare cases 
as low as 7200 A), where the deredshifted SDSS spectra do 
not extend up to this value. 

As opposed to the method used for the determination 
of the line intensities, for this population synthesis we 
used SSP spec t ra from the so-called " BC99 " SSP hbrary 
jW estera etall |2004| : ICuisinier et all I2OO6I : iLisker et al] 
200l|). It was produced using the GISSEL code^ imple ment- 
ing the Padova 2000 isochrones ('Girardi et al.' 2000") com- 
bined with the BaSeL 3.1 "Pado va 2000" stellar library 
(|Westera et al.ll2002l : IWesterall200ll ). The BaSeL 3.1 library 
was calibrated to reproduce the spectral shapes of stars of 
metallicities [Fe/H] from -2 to 0, so it is the ideal choice 
for full spectral fitting. By contrast, the "BC03" library 
described in the previous section would not have been a 
good choice, for the reasons mentioned there. As the GIS- 
SEL spectra do not include nebular continuum emission, 
w e added it to the sp e ctra, in the sarne way as described 
inlWestera et all (|2004l ): ICuisinier et al.l (|2006l ): lLisker et"al] 
(|2006l ). 

In lWestera et all l|2004h . we found that the stellar con- 
tent of HII galaxies is made up of at least three populations, 
a young one (up to 10 Myr old) that is responsible for the 
ionisation of the gas and, therefore, for the nebular emission, 
an intermediate one (from 20 to 500 Myr), and an old one 
(at least 1 Gyr old). We showed that these three popula- 
tions are necessary - and sufficient - to characterise an HI I 
galaxy, which is now the current view (|Hovos et al.ll2007l ). 
Therefore, in this work we modelled the actual population 
as being composed of an old, an intermediate, and a young 
stellar population. In a full spectral fit, we determined the 
masses, ages and metallicities of these partial populations. 
To obtain meaningful results, we chose to reduce the num- 
ber of fitting parameters (and thus, of degeneracies) to a 
minimum. 

As the young population should be more or less coeval 
with the gas, whose metallicity is the one attained now 
by these galaxies, we assumed the young and intermedi- 
ate populations to have the same chemical composition as 
the gas. As the gas metallicities we measured are in terms 
of [O/H], whereas our stellar population library was cal- 
ibrated in [Fe/H], we still had to translate the oxygen - 
into iron abundances, since the element ratios in HII galax- 
ies are far from solar. Oxygen is released by stellar winds 
and by SNe I and II, whereas iron and other iron peak 
elements are only produced in SNela, which appear only 
around 1 Gyr after the formation of a stellar population. 
Therefore, the ratio between the abundances of these two 
elements in a galaxy depends on its entire star formation 
history. Adding to this the possibility of "selective mass 
loss", as is predicted by the models mentioned in the in- 
troduction, it would be illusory to expect a unique iron-to- 
oxygen ratio for the galaxies of our sample. However, as our 
sample is very homogeneous, we do expect there to exist a 
mean relation and only a moderate scatter around this rela- 
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tion. Unfortunately, transformations between the two metal- 
licity indicators for HII galaxies are rare in the literature. 
In earlier abundance st udies of this type of gala xies, such 
as iKniazev et all (|2004l ') and llzotov et ahl (|2006l '). the iron 
abundance is not determined, as they are also based on emis- 
sion lines in SDSS spectra. On the other hand, most theoret- 
ical transformations, like the ones determin ed for the Milky 
Way b ulge and the solar neighbourhood bv lMatteucci et all 
1 19991 ). are not valid for HII galaxies. Finally, we derived the 
[O/F e] to [Fe/H] re l ation from the one for irregular galaxies 
from ICalura et all (|2009l 'l. who calculated the element ra- 
tios of galaxies of different morphological types using chem- 
ical evolution models. According to their fig. 13, the [O/Fe] 
to [Fe/H] relation for irregular galaxies is around [O/Fe] = 
-0.25 X [Fe/H] -0.425 ±0.15 independent of redshift, which 
translates to [Fe/H] = 1.333 x [0/H] -I- 0.5666 ± 0.2. Hence, 
we fixed the metaUicities of the young and intermediate pop- 
ulations, [Fe/H]y and [Fe/H]i, to 



[Fe/H] = [Fe/H]; = 1.333 x [0/H]p + 0.5666. 



(3) 



For the galaxies of our sample, most of which have [0/H]p 
values in the range from -1 to -0.5, we obtain [O/Fe] ratios 
between -0.2 and -0.4. As opposed to the metallicity of the 
young and intermediate populations, the metallicity of the 
old population, [Fe/H]o, is a free parameter of our fitting 
procedure, limited to the metaUicities of the spectra of the 
"BC99" SSP library, -2.252 -1.65, -1.25, -0.65, -0.35, 0.027, 
and 0.225. Contrary to our approach in previous works, we 
did not force the metaUicities of the old populations to be 
lower than the ones of the younger stars, since one of the 
aims of this work is to verify if they really are. 

As the decomposition of galaxy spectra into SSP spec- 
tra is known to present various degeneracies, we verified if 
our fitting procedure is able to recover the population pa- 
rameters of synthetic composite spectra made up of three 
SSPs with noise added. We found that the masses, ages, 
and metaUicities of the input populations were properly re- 
covered for the signal-to-noise ratios of both our HII galaxy 
sample and our test sample (^8 resp. 5 in the g band), and 
conclude that the number of free parameters of our proce- 
dure is adequate and that the results presented in the fol- 
lowing can be trusted to be meaningful. 

In order to determine the total stellar masses of our 
galaxies, we had to know their luminosity distances _Dl, and 
their (stellar) mass-to-light ratios. We calculated the former 
from the redshifts given by the SDSS, using a cosmology 



of Q,n 



0.3, Qa = 0.7, and Ho 



70 



km 
s-Mpc 



The steUar 



mass-to-light ratios result from our best fits. 



4 RESULTS 

The strengths of the emission lines used for this analysis, 
as well as the internal reddening values E{B — V), the gas 
metaUicities [0/H]p and the total stellar masses A/tot, can 
be found in Table [XT] of the electronic version of this paper. 

Fig. 0] shows the stellar mass and gas metallicity dis- 
tributions of our sample (in black) and the test sample (in 
grey). The average of the logarithms of the stellar masses 
(in Mq) amounts to 9.4 and the average gas metallicity to 
-0.79, whereas those values are 9.8 resp. -0.78 for the test 
sample. 




[0/H]p 

Figure 4. Stellar mass (top panel) - and gas metallicity (bottom 
panel) distributions of the galaxies of our sample (in black). The 
mass - and gas metallicity distributions of the test sample are 
shown in grey (scaled to the same total number of galaxies). 



As expected, the galaxies of the HII galaxy sample have, on 
average, lower masses than the ones of the test sample. Since 
our sample contains very few galaxies with stellar masses 
below 2 • 10* Mq or above 2 • 10^° Mq (The mass bins at 
5 ■ 10*^ Mq, 1 ■ 10^ Mq, 2 ■ lO'' Mq, 5 ■ lO'' Mq, 1 ■ 10** Mq; 
5 ■ 10^° Mq, 1 • 10" Mq, 2 • 10" Mq, 5 • 10" Mq, and 
1 • 10^^ Mq contain 1, 0, 4, 2, 3; 1, 2, 0, and galaxies, 
resp.), statistical studies of galaxies in these mass ranges 
suffer from low number statistics. In the following we will 
only consider the range from 2 ■ 10* Mq to 2 • 10^'' Mq when 
studying galaxy properties as a function of mass. In Figs. [7] 
and [8] this mass range is delimited by dotted lines. For the 
test sample, the mass range from 2 ■ 10* Mq to 1 ■ 10^^ Mq 
could be used for such studies, as each mass bin contains 
over 150 galaxies. 

The average gas metaUicities are nearly the same for both 
samples, but the distribution is narrower for the HII galax- 
ies, indicating that this is a more homegeneous sample. 

Fig. [5] shows the metallicity distributions of the par- 
tial populations. As we fixed [Fe/H]^ = [Fe/H]j = 1.333 x 
[0/H]p -(-0.5666 in the fitting procedure, we show in the up- 
per panel the metallicity of the joint young and intermediate 
age population. The average metallicity of the joint young 
and intermediate populations is -0.49, and of the old one 
-0.87 (test sample: -0.49, and -0.84). Here too, the distri- 
butions are narrower for the more homegeneous HII galax- 
ies sample. Nevertheless, the spread in [Fe/H]^ is wide. As 
[Fe/H]^ is a measure of the average metallicity of all pop- 
ulations older than ~ 1 Gyr, this wide spread refiects the 
various stages of (metallicity) evolution of the galaxies of 
our sample around one Gyr ago. 
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Figure 5. Metallicity distributions of the partial populations in 
our galaxy sample. Top panel: the joint young and intermediate 
age populations, bottom panel: old populations. The metallicity 
distributions of the partial populations in the test sample are 
shown in grey (scaled to the same total number of galaxies). 



Figure 6. Gas metallicities as a function of the metallicities of 
the old populations. The crosses represent the average of the gas 
metallicities for each old population metallicity value, and the 
error bars show the standard deviations. 



5 DISCUSSION 

In Fig. |6l we show the gas metallicity in function of the 
metallicity of the old population, both for the HII galaxies 
and for the test sample. In either sample, not much of a 
relation can be seen, any possible trend being significantly 
smaller than the scatter. Apparently, the stage of a galaxy's 
evolution around one Gyr ago is not necessarily reflected in 
the present day one. The chemical evolution and the present 
day (gas) metallicities of galaxies suffer stronger influences 
from other factors, such as the galaxy mass. 

Fig. [7] shows, how the present gas metallicity [0/H]p 
depends on the total stellar galaxy mass. In the test sam- 
ple, [0/H]p clearly increases with mass in the range from 
2-10* Mq to 1- 10"*^"*^ M(7), a s has been found in previous stud- 



ies (Tremonti et al 



I2OO5I : iPanter et al 



.20041 : lAsari et all l2007l : iThomas et al] 
20081 . and others). However, our mass- 



metallicity relation lies lower than the ones found in these 
works, as recent studies (|Ellison et al.|[200^ : lMannucci et al.l 
l2010h suggest to be the case for samples of star forming 
galaxies, such as our test sample. 

The HII galaxy sample represents galaxies of a more 
homogeneous class than the test sample. Therefore, in 
this sample, the mass-metallicity relation is much less pro- 
nounced, if existing at all. In fact, the average gas metallic- 
ity stays constant between -0.7 and -0.8 in t he full galaxy 
mass range from 2 • 10* M© to 2 • 10^° Mq. IPanter et al] 
(|2008l ) find almost the same relation for their sub-sample 
of galaxies that are (flux-) dominated by young populations 
(^ 0.5 Gyr), as can be seen in their fig. 10. 
In a "closed box" context, this would mean that the gas 
fraction that has been transformed into stars so far is in 
the same range for all galaxies of the HII galaxy sample. 
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Figure 7. Top panel: the gas metallicity [0/H]p as a function of 
the total stellar mass of the galaxies of our sample. The dots show 
the individual galaxies, whereas the crosses, boxes, and bars show 
the statistical properties after subdividing the sample into mass 
bins. For each value of Mtot, the crosses represent the average, the 
boxes represent the standard deviation, and the vertical bars rep- 
resent the 90th percentile range. The dotted vertical lines detach 
the "trustworthy" range from the nearly empty mass bins. Bot- 
tom panel: the individual galaxies, average, standard deviation, 
and 90th percentile range of [0/H]p for the test sample. 
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Figure 8. Top panel: the metal enrichment, [Fe/HJ^^j — [Fe/H]^, 
as a function of the total mass of the galaxies of our sample. For 
each value of Mtot, the crosses represent the average, the boxes 
represent the standard deviation, and the vertical bars represent 
the 90th percentile range. The dotted vertical lines detach the 
"trustworthy" range from the nearly empty mass bins. Bottom 
panel: the average, standard deviation, and 90th percentile range 
of [Fe/H]y/; - [Fe/H]^ for the test sample. 



independent of their masses. However, HII galaxies are not 
"closed boxes" , and our results should be interpreted in the 
light of models including interaction processes with the en- 
vironment, like the ones presented in the introduction. 
Most of these models predict that HII galaxies lose a large 
fraction of their heavy elements through galactic winds, 
whose escape efficiencies depend upon the galaxy mass. 
These models are at odds with our result for two reasons: 
They result in a dependence of the metallicity on the galaxy 
mass, and their predicted present day metallicities are lower 

than the ones we measu re. 

The iRecchi et all (|2006l ): iRecchi fc HensleJ (|2007l ) models 
with constantly infalling clouds ap pear more promisin g, es- 
pecially th e models NGC-4-1- B C of IRecchi et all (120061 ) and 
NCSM of IRecchi fc Hensleij (120071 ). The clouds partially 
retain the metal-rich gas, resulting in higher present-day 
metallicities. Depending on the star formation history and 
on the mass of the infalling clouds, galaxies like the ones 
described in these models can reach the metallicities we 
measure. On top of that, this scenario provides a possible 
explanation, why we find metallicities in the same ranges 
for different stellar masses. Not only the SFR and, there- 
fore, the amount of produced metals and galactic wind in- 
creases with increasing galaxy mass, but also the amount 
of infalling low- metallicity clouds and of retained outflowing 
high-metallicity gas. Of course, only detailed model calcu- 
lations can tell, if the mass dependencies of the different 
metallicity-increasing and -decreasing processes really add 
up in such a way, that the final metallicity becomes mass 
independent. The fact, that a constant mass-metallicity re- 



lation is observed only in our HII galaxy sample, and not, 
for example, in our test sample, or in any of the studies 
mentioned at the beginning of this section, shows that only 
in galaxies with extreme star formation, and thus extreme 
stellar winds, the metal outflow manages to balance the re- 
tention in this way. However, the models including infalling 
clouds seem to be on the right track. 

Another noteworthy point of Fig. [7] is that, though the av- 
erage gas metallicity [0/H]p does not vary with the galaxy 
mass, its dispersion does increase with galactic mass. This 
indicates that high mass galaxies have a more complex chem- 
ical, and certainly assembly, history. 



Apart from the present-day metallicity, our analysis also 
allows us to quantify the chemical evolution of our sample 
galaxies since the formation of the old population. As the 
metallicity of the young (and intermediate) population rep- 
resents the galaxy's present-day metallicity, and the one of 
the old population is a measure of the average metallicity of 
all stars older than ~ 1 Gyr, the quantity [Fe/HJ^yj — [Fe/H]^, 
which we shall call the metal enrichment, is a good measure 
of the change in metallicity during the last few Gyr. Fig. [8] 
shows the metal enrichment in function of the stellar mass 
for the galaxies in both our HII galaxies sample and the 
test sample. For our HII galaxies, the enrichment is fairly 
constant of the order of 0.3 up to 5 • 10® Mq, and almost 
0.5 dex higher for the last two mass bins within the mass 
range, where reliable statistics are possible, 1 ■ 10^" Mq and 
2 • 10^" Mq . This suggests that the enrichment mechanism 
are independent of mass up to 5 • 10® Mq, which is com- 
patible with the formation scenario favoured in the previous 
paragraph. For galaxy masses above 5 • 10® Mq, the en- 
richment is significantly higher, which means that at high 
masses, the metal-rich gas could not escape the galaxy's po - 
tential well, as predicted in the iMacLow fc Ferraral (Il999l ) 
models. Here too, the dispersion increases with galaxy mass, 
confirming more complex chemical and assembly histories 
for high mass galaxies. 

An interesting point is that we find a negative metal enrich- 
ment, thus a metal reduction, for 164 galaxies, that is, 23 per 
cent of our sample. Two mechanisms are known to reduce 
the (gas) metallicity of a galaxy: Selective outflow of metal- 
enriched gas and infall of low-metallicity gas clouds, maybe 
even of primordial composition. Both processes might have 
contributed to the metal reduction of our galaxies, but in 
the light of our previous results, we expect the second one 
to be more important on long time-scales. 



It is intriguing that in the test sample, too, we find 
the metal enrichment to be mass independent, even though 
the final metallicity increases with galaxy mass. [Fe/HJ^^j — 
[Fe/H]^ amounts to around 0.36 in the full range from 
2 • 10** Mq to 1 • 10^^ Mq. Even the dispersion is fairly con- 
stant over this mass range. Apparently, enrichment mecha- 
nisms are similar for galaxies of different masses even in a 
heterogeneous sample like our test sample. Around 44 per 
cent of the test sample galaxies, 33178 in total, show nega- 
tive metal enrichment, which proves that processes like se- 
lective outflow or low metallicity gas infall are also at play 
in these galaxies. 
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6 CONCLUSIONS 

We performed a gas - and stellar population analysis of HII 
galaxies from the Sloan Digital Sky Survey Data Release 7, 
using their full spectra. We selected a sample which contains 
712 HII galaxy spectra, a relatively low number due to our 
high requirements towards the quality of the spectra and to 
the criterion that the spectra contain all the spectral lines 
necessary for our study. We also selected a test sample of 
74989 galaxies, using less restrictive requirements, the only 
criterion being that they show some line emission. There- 
with, it is a sample of star forming galaxies. 
We derived independent metallicities for the young popula- 
tions present in our sample galaxies from the gas emission 
lines, and for the old populations from full spectra fitting, 
that is, from the continuum and absorption features from 
stars. 

Our test sample follows the well-known mass-metallicity 
relation, e.g. the fact that the (gas) metallicity increases with 
(stellar) galaxy mass, indicating that low mass galaxies are 
chemically less evolved than high mass galaxies. We do not 
find any systematical tendency of the difference between the 
metallicities of the young and the old stellar populations, a 
quantity we call metallicity enrichment, with galactic mass, 
indicating that, on average, recycling mechanisms of the in- 
terstellar medium should be the same for low - and high 
mass galaxies. 

In our HII galaxy sample, on the other hand, we find 
no mass dependence of the present day metallicity. For 
the metallicity enrichment, we find no mass dependence for 
masses up to 5 • 10^ Mq, whereas above this mass, the en- 
richment is significantly higher. 

We interpret these findings in the light of recent hy- 
drodynamical evolutionary models of dwarf irregular galax- 
ies, of which HII galaxies are a sub-class. We favour mod- 
els featuring constantly infalling low-metallicity clouds able 
to retain part of the high metalli city galactic winds, such 
as the mode ls NGC-4-I-BC fro r ti iRecchi et al.l (|2006l ) and 
NCSM from iRecchi fc Hensled l|2007h . since models with- 
out infalling clouds fail to predict the metallicities and the 
non-dependence on galaxy mass of our sample galaxies. For 
galaxies with masses above 5 • 10^ Mq, most, if not all, of 
the metal-rich gas is reta ined by the galaxy's gravitational 
potential, as predicted bv lMacLow fc Ferraral l|l999r ). which 
leads to the observed break-down of the constant mass- 
enrichment relation. 

However, the dispersion of both the present day metal- 
licity and the metallicity enrichment does increase with 
galactic mass, indicating more complex chemical and assem- 
bly histories for high mass HII galaxies. 
In 23 per cent of our HII galaxies, the metallicity has de- 
creased in the last few Gyr. We attribute the metallicity 
decrease mainly to the dilution of the galactic gas by the 
infalling clouds, which is compatible with the models we 
favour. 
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APPENDIX A: POPULATION PARAMETERS 
OF INDIVIDUAL SPECTRA 

This paper has been typeset from a TjnjX/ M^gX file prepared 
by the author. 
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Table Al. Line strengths, internal reddening, gas metallicities, and total stellar masses of individual spectra (first 10 lines only, for the 
complete sample, see the electronic edition). 



SDSS spectroscopic ID 


/(H/3)* 


H«** 


[Oil]" 




[OIII]" 


[OIII]" 


[Nil]" 


[Nil]" 


E{B-V) 


[0/H]p 


Mtot 


MJD-plate-fiber 


4861 A 


6563 A 


3727 A 


3727 A 


4959 A 


5007 A 


6548 A 


6584 A 






[109 A/q] 


51637-0306-583 


5.102 


3.010 


5.848 


2.948 


0.728 


2.088 


0.148 


0.472 


0.088 


-1.022 


0.63 


51658-0282-047 


10.075 


3.103 


5.644 


2.806 


0.694 


2.049 


0.143 


0.457 


0.145 


-0.991 


1.93 


51662-0308-628 


10.264 


3.017 


4.800 


4.074 


0.669 


2.049 


0.158 


0.571 


0.093 


-0.860 


2.11 


51663-0307-268 


14.076 


3.052 


4.270 


3.005 


1.009 


3.078 


0.176 


0.574 


0.114 


-0.794 


3.27 


51671-0299-571 


8.794 


3.179 


6.656 


2.846 


0.615 


2.199 


0.203 


0.626 


0.189 


-1.142 


6.78 


51691-0350-439 


24.324 


3.105 


3.533 


3.102 


0.550 


1.696 


0.212 


0.644 


0.146 


-0.654 


1.98 


51692-0339-437 


11.115 


2.921 


2.707 


4.272 


1.238 


3.741 


0.074 


0.223 


0.033 


-0.599 


0.34 


51783-0395-570 


9.257 


2.924 


3.362 


4.199 


0.639 


1.975 


0.118 


0.374 


0.034 


-0.624 


0.98 


51812-0404-507 


9.680 


3.042 


4.195 


3.339 


0.641 


1.964 


0.144 


0.490 


0.108 


-0.763 


4.10 


51818-0383-266 


8.743 


2.944 


3.631 


3.398 


0.899 


2.789 


0.097 


0.303 


0.047 


-0.687 


0.87 



*/{H/9) given in units of 10"*^^ erg/s/cm^ 
**flux ratios given in /(A)/7(H/3) 
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